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ABSTRACT 

In  ultrasonic  NDE,  simulation  studies  can  play  an  important  role  in 
complimenting  experimental  validation  of  techniques  under  development.  The  utility  of 
such  simulations  depends,  in  part,  on  the  degree  to  which  the  simulated  defect  and  noise 
signals  are  representative  of  the  measured  signals.  In  this  paper,  we  describe  an  approach 
for  generating  simulated  acoustic  noise  with  a  spatial  correlation  coefficient  distribution 
and  maximum  extreme  value  (MEV)  distribution  which  matches  those  distributions  for 
measured  acoustic  noise.  The  procedure  for  generating  noise  signals  is  outlined  for  a  line 
scan  and  for  a  raster  scan.  The  basic  approach  forces  the  correlation  of  neighboring 
signals  to  the  desired  correlation  by  creating  each  signal  as  the  sum  of  appropriately 
scaled  neighboring  signals  plus  a  new  random  signal.  Eor  the  line  scan  where  each 
interior  position  has  only  two  neighbors,  this  process  is  done  sequentially  without 
iteration.  Eor  the  raster  scan  where  each  interior  point  has  four  nearest  neighbors, 
iteration  is  required  to  simultaneously  achieve  the  desired  correlations  with  row  and 
column  neighbors.  The  MEV  distribution  is  controlled  in  an  outer  iterative  loop  with  the 
shape  and  position  of  the  distribution  dictated  by  spectral  content  of  the  noise  signals  and 
by  controlling  the  signal  energy,  respectively.  Results  are  shown  which  demonstrate  the 
effectiveness  of  the  approach.  With  this  approach,  a  limited  number  of  measured  signals 
can  be  used  to  establish  the  correlation  coefficient  and  MEV  distributions  which  drive  the 
computer  generation  of  a  large  number  of  simulated  acoustic  noise  signals. 

Keywords:  ultrasonics,  acoustic  noise,  correlation 

PREPRINT 


I,  INTRODUCTION 


Simulation  studies  are  routinely  used  in  ultrasonic  nondestructive  evaluation 
(NDE)  to  complement  experimental  studies  during  the  development  of  new  inspection 
approaches.  Simulated  A-scans  typically  include  a  target  signal  (e.g.,  a  flaw  or  weld 
plane  signal)  plus  noise.  The  noise  of  interest  lies  between  the  front  and  back  surface 
reflections  in  an  A-scan  and  is  comprised  of  electronic  plus  acoustic  noise.  In  certain 
cases,  the  degree  of  correlation  between  A-scans  can  have  a  significant  influence  on 
detection.*^''"^^  These  correlations  can  be  quantified  in  terms  of  the  spatial  cross¬ 
correlation,  at  zero  lag,  between  gated  A-scans. Uncorrelated  noise  A-scans  are  easily 
generated  using  a  normal  random  number  generator  with  filtering  used  to  achieve  the 
desired  frequency  content.  Margetan  et  al.  have  gone  a  step  further  using  an  independent 
scatterer  model  to  generate  simulated  grain  noise  A-scans  with  correlations  between  A- 
scans  controlled  to  some  degree  through  use  of  an  experimentally  determined  spatial 
correlation  length  parameter.*^^^  Using  their  approach,  the  average,  maximum,  and 
standard  deviation  of  gated  peak-values  are  also  controlled.  In  the  current  paper,  we 
describe  an  approach  for  generating  simulated  acoustic  noise  with  a  correlation 
coefficient  distribution  and  maximum  extreme  value  (MEV)  distribution  (i.e,  the  gated 
peak-value  distribution)  which  matches  those  distributions  determined  for  measured 
acoustic  noise. 

The  motivation  for  this  project  finds  its  origin  in  research  associated  with  kissing 
bond  detection  for  inertial  welding  of  two  stainless  steel  pieces.  Kissing  bonds  can  be 
very  difficult  and  expensive  to  fabricate  in  a  controlled  fashion,  providing  motivation  for 
simulation  studies  which  utilize  simulated  noise  signals  with  realistic  correlation 
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coefficient  distributions.  The  basic  approach  being  developed  relies  on  the  correlation 
coefficients  between  adjacent  A-scans  to  detect  low  signal-to-noise  ratio  (SNR)  kissing 
bond  signals.  This  approach  finds  its  genesis  in  the  work  of  Nagy  and  Adler  on  this  same 
problem.*^^^  Without  going  into  the  analysis  details,  suffice  it  to  say  that  an  inspection 
approach  can  be  formulated  which  relies  heavily  on  the  comparison  between  correlation 
coefficient  distributions  associated  with  backscattered  signals  from  the  weld  being 
inspected  and  from  a  known  set  of  acceptable  welds. More  recently  the  correlation 
approach  has  been  extended  to  crack  detection  in  both  pulse/echo  and  pitch/catch. 

Computer  generation  of  correlated  random  variables  with  a  desired  mean 
correlation  coefficient  is  straightforward.  For  example,  the  random  sample  can  be 
generated  based  on  the  -1  random  sample  as  follows: 

x{i)  =  a{i)  + bx{i  - 1)  (1) 

where  b  is  an  adjustable  scale  factor,  a(z)  is  the  output  of  a  random  number  generator, 
and  the  x 's  are  the  computer  generated  samples  of  the  random  variable,  x.  The  single 
scale  factor  can  be  adjusted  to  yield  the  desired  mean  correlation  coefficient  value,  but,  in 
general,  the  scale  factor  cannot  be  adjusted  such  that  desired  correlation  coefficient  mean 
and  distribution  (width  and  shape)  are  achieved.  The  natural  relationship  between  the 
mean  and  shape  dictate  that  as  b  is  increased,  correlation  values  are  forced  closer  to  the 
limiting  value  of  1.0,  the  distribution  breadth  decreases,  and  the  distribution  becomes 
increasingly  skewed  with  a  fat  lower  tail.  The  primary  task  addressed  in  this  paper  is  to 
extend  the  approach  represented  in  Eq.  (1)  so  that  vectors  of  random  numbers  (simulated 
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A-scans)  can  be  generated  that  show  the  desired  distribution  of  eorrelation  eoeffieient 
values.  These  simulated  signals  will  also  be  foreed  to  the  desired  MEV  distribution. 

Simulation  of  data  with  speeified  eorrelation  strueture  arises  in  many  fields.  In 
the  physieal  seienees,  seientists  often  use  Shewhart  Control  Charts  to  monitor  a  proeess. 
In  this  setting,  data  are  a  time-indexed  series  of  averages  or  eounts  taken  at  regularly 
spaeed  time  intervals.  Padgett,  Thombs  and  Padgett  present  a  method  to  generate  sueh 
one-dimensional  data  with  speeified  mean  and  varianee  strueture  so  that  the  performanee 
of  sueh  eharts  ean  be  studied  via  simulation.  In  this  paper,  the  generated  data  are  two 
dimensional,  with  emphasis  on  both  the  eorrelation  within  and  between  series,  as  well 
and  the  maximum  extreme  value  distribution. 

Two  dimensional  spatial  time  series  data  with  eorrelation  strueture  is  eommon  in 
many  other  areas,  ineluding  meteorology,  hydrology  and  eeologieal  and  environmental 
studies.  In  the  area  of  eeology  and  wildlife  studies,  observations  sueh  as  animal  eounts 
are  typieally  observed  at  irregularly-spaeed  loeations.  Neighboring  observations  are 
likely  to  be  eorrelated,  and  the  paper  by  Brooker  represents  one  of  the  first  attempts  to 
generate  sueh  data. 

In  spatial  statisties,  point  patterns  are  observed  on  variables  sueh  as  temperature 
and  rainfall,  so  that  the  data  often  inelude  a  (third)  eomponent,  time.  See  Cressie  for 
more  information  on  the  statistieal  aspeets  of  fitting  models  to  sueh  data.  A  reeent 
eontribution  by  Kyriakidis  et  al.^'^^  proposes  an  algorithm  for  generating  spatio-temporal 
preeipitation  data,  with  emphasis  on  preserving  the  distribution  of  the  original  data  set. 
Data  are  both  spaee-indexed  (e.g.,  longitude  and  latitude)  and  time  indexed. 
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The  methodology  presented  in  the  eurrent  paper  is  distinet  from  these  related 
approaehes  for  generating  dependent  data  in  that  both  the  eorrelation  and  the  MEV 
distribution  are  eontrolled.  The  paper  proeeeds  by  first  establishing  representative 
eorrelation  eoeffieient  and  MEV  distributions  assoeiated  with  measured  baekseattered 
noise.  The  methodology  for  generation  of  simulated  acoustic  noise  with  the  desired 
correlation  coefficient  and  MEV  distributions  is  then  described  for  a  line  scan  and  for  a 
raster  scan.  Results  are  presented  which  validate  the  approach  for  the  experimentally 
established  correlation  coefficient  and  MEV  distributions.  The  paper  closes  with  a  brief 
summary  section. 

II,  MEASUREMENT  PROCEDURE 

Baekseattered  grain  noise  signals  were  measured  and  used  as  a  basis  for 
calculating  associated  correlation  coefficient  and  MEV  distributions.  The  ultrasonic 
measurement  system  used  in  these  measurements  consists  of  a  water  tank  filled  with 
degassed  tap  water  at  approximately  19°C,  a  three  dimensional  scanning  bridge  that  holds 
a  transducer,  a  pulser-receiver  unit,  and  a  12  bit  data  acquisition  card  with  a  sample  rate 
of  100  MS/s.  A  dedicated  PC  collects  data  from  the  acquisition  card  and  controls  the 
motor  controller  that  moves  the  scanning  bridge.  A  separate  PC  is  employed  for  data 
analysis. 

The  transducer  used  to  make  measurements  was  a  focused  14”  transducer  with  a 
10  MHz  center  frequency  and  a  4”  focal  length.  The  settings  on  the  pulser-receiver  and 
the  data  acquisition  card  were  typically  set  at  values  such  that  the  front  and  back  surface 
reflections  were  blown  off  the  screen  in  order  to  enable  proper  digitization  of  the  grain 
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noise.  In  each  measurement  position,  64  signals  were  taken  and  averaged  together  in 
order  to  reduce  electronic  noise. 

The  sample  used  in  all  measurements  was  a  stainless  steel  plate  with  the 
dimensions  10.1  x  5  x  1.9  cm.  A  leveling  plate  was  used  to  ensure  that  the  specimen  was 
aligned  with  the  transducer’s  scan  plane.  The  transducer  was  normalized  in  relation  to 
the  front  face  of  the  sample.  Data  was  taken  at  a  single  water  path  such  that  the  focal 
point  of  the  transducer  would  be  approximately  at  the  sample  mid-plane.  Signals  were 
measured  on  a  3.2  x  2.4  cm  grid  with  0.5  mm  between  measurement  positions. 

After  the  data  was  collected,  pre-processing  was  done  on  the  raw  signals.  First, 
the  front  surface  reflections  of  each  signal  were  aligned  with  one  another.  Then  all  of  the 
signals  were  averaged  together  in  order  to  identify  any  non-random  component 
associated  with  the  front  surface  reflections.  This  mean  signal  was  subtracted  from  each 
of  the  individual  signals  so  that  spatial  correlations  between  adjacent  noise  signals  could 
be  calculated  with  minimal  influence  from  front  surface  reflection  ringing. Finally, 
each  signal  was  gated  to  extract  a  time  window  of  51  points  (0.50  ps)  with  the  gate 
positioned  to  correspond  closely  with  the  location  in  the  sample  where  the  transducer  was 
focused. 

Ill,  STATISTICAL  ANALYSIS  OF  MEASURED  NOISE 
A,  Maximum  Extreme  Value  Distributions 

Vectors  of  data  representing  maximum  extreme  values  (maximum  absolute 
values)  for  the  measured  noise  were  extracted  from  the  gated  signals.  Figure  1  shows  2 
histograms  based  on  this  MEV  data.  The  lower  histogram  is  for  all  measured  signals, 
that  is,  for  a  measurement  position  spacing  of  0.5  mm.  The  upper  histogram  is  for  every 
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3rd  signal,  corresponding  to  a  measurement  spaeing  of  1.5  mm.  The  histograms  are 
fundamentally  the  same  with  the  number  of  observations  being  the  only  signifieant 
differenee.  An  extreme  value  distribution  based  on  the  absolute  value  of  the  differenee 
between  min  and  max  values  eould  have  been  utilized  instead.*^^^ 

B,  Spatial  Correlation  Coefficient  Distributions 

The  eorrelation  of  interest  is  the  spatial  eross-eorrelation  ealeulated  at  zero  lag 
between  gated  A-seans  measured  at  adjaeent  measurement  positions.  For  diseussion 
purposes,  eonsider  an  A  x  M  sean  with  the  A-sean  (eaeh  T  points  long)  written  into  the 
matrix  x  =  x{i,j,t)  i  =  \,N  j  =  l,M  t  =  0,7  - 1  (see  Fig.  2).  Calculation  of  the  eorrelation 
eoeffieient  between  A-seans  is  then  given  by  the  following  equation  where 
p  ~  /’(h  y)  h  ’  5  )  • 


Iffy  (/,  7)][^0'  +  4  A'  +  ’  0  -  O'  +  4  A'  +  )] 


1  2 

1  2 

1  ^ + 4  ’  j  ^  4  ’  0  ~  0  ^  4  ’  j  ^  4  )J 

V=ti  \ 

t=ti 

In  Eq.  (2),  the  summation  range  in  the  time-domain  defines  the  portion  of  the  signal  (the 
time  window  or  gate)  of  interest,  m  is  the  mean  value  ealeulated  over  the  gate,  and  5  is 
a  spatial  shift  parameter.  Throughout  the  paper,  t  is  used  as  a  diserete  index  referring  to 
the  temporal  direetion.  With  4=1  4=0’  eorrelations  are  ealeulated  between 

adjaeent  A-sean,  that  is,  A-seans  measured  at  the  and  y'^+l  positions  in  the  i‘^  sean 
row.  Similarly,  with  4=0  4=1’  column  eorrelations  ean  be  ealeulated.  Correlations 
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for  all  possible  adjacent  signal  combinations  in  an  x  M  raster  scan  can  be  established 
using  a  computation  loop  over  i  and  j  with  the  spatial  shift  applied  sequentially  to  i  and  j. 

In  order  to  define  desired  distributions  of  spatial  correlation  coefficients  for 
simulated  noise,  conditional  correlation  coefficient  distributions  were  established  based 
on  the  pre-processed  measured  noise  signals.  In  this  section,  we  begin  by  considering  the 
overall  distribution  of  correlation  coefficients.  Conditional  distributions  are  addressed 
below.  The  first  step  was  to  calculate  the  correlation  coefficient  (Eq.  2)  between  each 
signal  and  its  row  and  column  neighbors  (see  Fig.  2)  and  write  these  values  into  a  single 
vector.  This  calculation  was  done  separately  for  three  measurement  spacings;  0.5  mm 
spacing  (every  signal),  1  mm  spacing  (every  2nd  signal),  and  2  mm  spacing  (every  4* 
signal).  The  correlation  coefficient  histograms  for  these  measurement  spacings  are 
shown  in  Fig.  3.  The  influence  of  measurement  spacing  on  both  the  mean  and  the  shape 
of  the  correlation  coefficient  distribution  is  apparent.  Conditional  correlation  coefficient 
distributions  show  the  same  basic  characteristics  as  these  overall  distributions. 

Probability  density  function  fits  to  the  distributions  were  considered.  The  normal 
and  gamma  distributions  were  found  to  be  the  two  common  distributions  that  yielded  the 
best  fit  to  the  correlation  distributions.  The  normal  distribution  works  best  when  the 
average  correlation  value  is  near  zero  and  the  distribution  is  symmetric;  however,  as  is 
apparent  from  the  figure,  the  distribution  tends  to  show  a  suppressed  peak  and  fatter  tails 
than  the  normal  distribution.  The  gamma  distribution  works  well  for  skewed 
distributions  with  higher  average  correlations.  In  some  cases,  a  standard  distribution 
which  shows  reasonable  fit  to  the  histogram  cannot  be  found;  however,  an  interpolation 
of  the  histogram  can  be  used  to  approximate  the  actual  distribution  of  the  correlation 
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coefficient  values.  The  interpolation  approach  utilizes  the  Matlab  routine  randsample  to 
randomly  choose  a  bin  and  then  randomly  select  a  value  from  each  bin  based  on  a  linear 
probability  density  function  defined  between  the  edges  of  each  bin.*^^^^  Superimposed  on 
each  histogram  in  Figure  3  is  a  normal  or  gamma  distribution  fit  along  with  the 
interpolation  fit. 

In  general,  correlation  coefficient  values  are  not  randomly  distributed  in  space. 
That  is,  the  correlation  coefficient  between  a  given  pair  of  signals  depends  on  the 
correlation  coefficients  between  surrounding  pairs  of  signals.  As  an  example.  Fig.  4 
shows  the  variation  in  correlation  values  along  a  line  scan  with  relatively  low  (upper 
graph)  and  high  (lower  graph)  measurement  spacing.  As  demonstrated  in  the  figure,  the 
dependence  of  correlation  values  on  adjacent  values  decreases  with  increasing 
measurement  spacing.  A  graphical  example  of  a  conditional  histogram  is  shown  by  the 
stem-plot  in  Fig.  5.  The  overall  distribution  is  given  by  stems  terminating  in  light  circles, 
and  the  conditional  distribution,  assuming  an  adjacent  correlation  value  of  0.58,  is  given 
by  the  stems  terminating  in  dark  circles. 

IV.  GENERATION  OF  SPATIALLY  CORRELATED  NOISE 
A.  Generation  of  Spatially  Uncorrelated  Acoustic  Noise 

A  number  of  equally  effective  approaches  could  be  taken  to  generate  uncorrelated 
acoustic  noise.  As  depicted  in  Fig.  6,  the  basic  steps  in  the  process  as  implemented  here 
are  as  follows:  1)  using  a  standard  normal  random  number  generator,  create  a  time- 
domain  white  noise  signal,  T  points  long,  for  each  (simulated)  measurement  position;  2) 
Fourier  transform  each  signal  to  the  frequency  domain;  3)  filter  each  signal  using  a  filter 
(with  unit  energy)  that  will  force  the  resultant  noise  to  have  the  same  average  power 
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spectrum  as  the  measured  noise;  and  4)  inverse  Fourier  transform  each  signal  back  to  the 
time-domain  to  yield  A-scans  which  are  spatially  uncorrelated  with  appropriate  frequency 
content.  Note  that  the  correlated  noise  creation  step  described  below  involves  sums  of 
scaled  A-scans.  This  process  changes  the  average  power  spectrum  from  the  desired 
spectrum  associated  with  these  initial  spatially  uncorrelated  signals.  An  iterative 
correction  approach  will  be  outlined  which  brings  the  average  power  spectrum  back  to 
the  desired  shape  while  still  achieving  the  desired  correlation  distribution  was 
implemented. 

B,  One  Dimensional  Generation  of  Spatially  Correlated  Noise 

We  begin  by  establishing  the  approach  for  a  line  scan  and  then  expand  the 
procedure  for  application  to  an  xy  raster  scan  in  the  next  section.  Starting  at  the  noise 
measurement  stage,  assume  that  a  line  scan  is  performed,  acoustic  noise  signals  are 
measured  and  pre-processed,  the  average  power  spectrum  is  estimated,  and  the  overall 
and  conditional  correlation  coefficient  distributions  are  established.  The  line  scan  of 
measured  signals  are  written  in  a  two-dimensional  matrix,  denoted  ,  with  average 
power  spectral  density  function  estimate  and  associated  magnitude  spectrum  given  by 
and  |X„(/)|,  respectively. 

For  discussion  purposes,  assume  that  an  N  position  line  scan  is  to  be  simulated. 
N correlation  coefficients  are  randomly  generated  with  the  first  value  coming  from 
the  overall  distribution  and  subsequent  values  drawn  from  conditional  correlation 
coefficient  distributions.  These  A-1  correlation  coefficients,  denoted  p{i)  i  =  2,N, 
will  dictate  the  correlations  between  simulated  noise  signals. 
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Clarification  of  the  overall  and  eonditional  eorrelation  eoefficient  distributions  is 


in  order.  The  overall  distribution  is  the  probability  density  funetion  assoeiated  with  the 
probability  /’(/?(/))  ,  that  is,  for  the  entire  set  of  eorrelation  coeffieients  given  that  these 
correlation  eoeffieients  are  not  independent.  Before  addressing  conditional  distributions, 
recall  the  notation:  p(i)  gives  the  eorrelation  between  x{i,t)  and  x(i-\,t)\  >o(t-l) 
gives  the  eorrelation  between  x{i-\,t)  and  x(z-2,t);  and  +  gives  the  eorrelation 
between  x{i  +  \,t)  and  x{i,t).  Considering  only  nearest  neighbors,  the  conditional 
probability  density  function  for  /?(/)  is  associated  with  the  probability: 
P{p{i)\  p{i -\)  and  /?(z  +  l))  .  In  practiee,  during  the  sequential  generation  of  eorrelated 
signals,  the  eorrelation  eoeffieient  p{i)  is  drawn  from  the  distribution  assoeiated  with 
sinee  p{i  + 1)  does  not  yet  exist. 

The  proeess  of  ereating  eorrelated  noise  signals  involves  several  steps.  We  begin 
by  deseribing  the  proeedure  used  to  foree  the  desired  correlation  between  to  two  given 
signals.  This  approaeh  is  then  ineorporated  into  an  iterative  proeedure  used  to 
simultaneously  mateh  the  desired  eorrelation  distribution  and  average  power  speetrum. 

Figure  7  gives  a  eartoon  representation  of  the  notation  and  some  of  the  steps 
involved  in  ereating  a  simulated  line  sean  of  spatially  correlated  signals.  Using  the 
proeedure  deseribed  in  the  previous  seetion,  the  first  step  is  to  ereate  a  set  of  N 
uneorrelated  aeoustie  noise  signals,  represented  by  the  matrix 
a  =  a(z,t)  i  =  \,N  t  =  0,7-1.  The  output  matrix,  i.e.,  the  set  of  eorrelated  signals,  is 
denoted  x  =  x(i,t)  i  =  \,N  t  =  0,7 -1 .  The  next  step  is  to  set  x(\,t)  equal  to  that 

is,  x(l,  t)  =  a(l,  t)  t  =  0, 7  - 1 .  Throughout  the  remainder  of  the  paper,  all  operations  will 
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be  implicitly  carried  out  over  the  range  of  t,  for  example, 

x(l,t)  =  a(l,t)  ^  x(l,t)  =  a(l,t)  t  =  0,7-1.  The  second  signal,  is  then  to  be 

determined  so  that  the  correlation  between  x{l,t)  and  x(l,t)  is  equal  to  the  desired 

correlation,  >5(2).  Consistent  with  approach  of  Eq.  (1),  the  method  used  to  calculate 

x(2,t)  is  to  add  b{l)x{\,t)  to  the  uncorrelated  signal  a{l,t):  x{l,t)  =  a{\,t)  +  b{l)x{\,t) . 
The  key  is  to  find  the  value  of  5(2)  which  forces  the  correlation  between  x(\,t)  and 
x{2,t)  to  be  equal  to >5(2).  This  process  will  be  repeated  for  each  signal,  that  is: 

x{i,  t)  =  a{i,  t)  +  b{i)x{i  - 1,  t)  i  =  2,N  (3) 

Note  that  x{i,t)  will  be  appropriately  correlated  with  both  of  its  neighbors,  x(i-\,t)  and 
x{i  +  \,t),  since  the  process  forces  correlations  between  x{i,t)  and  x{i-\,t)  and  between 
x{i  +  \,t)  and  x{i,t) . 

The  process  for  establishing  the  scale  factor,  5(z),  can  be  described  as  follows. 
We  begin  by  using  the  correlation  coefficient  between  a{i,t)  and  x(z-l,t)  to  calculate 
the  inherent  similarity  between  the  two  signals  that  will  be  used  to  create  the  output 
signal. 


T-\ 


^  [a{i,  t)  -  Ma  {i)){x{i  - 1, 0  ■ 

t=0 


T-l 

^(a(z,0 

t=0 


-mji, 


T-l 

^(x0'-l,0 

t=0 


i  =  2,N 


:\\2 


(4) 
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In  order  to  simplify  the  notation,  the  mean  subtraetion  step  will  be  implieit  in  future 
correlation  equations.  For  ease  of  calculations,  the  following  substitutions  will  be  used. 


T-l 

Ti  = 

T-\ 

^=0 

T-l 

^=0 


(5) 


The  goal  is  to  force  the  correlation  between  x(i,t)  and  x(i-l,t)  to  be  equal  to  p(i) .  This 
correlation  can  be  calculated  as  follows. 


T-l 

Y,x{i,t)x{i-\,t) 


t^O 


i 


T-l 

J^x(i,tf 

t=0 


T-l 

J^x(i-i,ty 

(=0 


(6) 


To  relate  this  correlation  to  the  scale  factor,  b(i),  we  proceed  with  manipulations  of  Eq. 
(3).  First,  both  sides  of  Eq.  (3)  are  multiplied  by  x(i  -  \,t)  to  give  the  top  equality  in  Eq. 
(7).  Second,  both  sides  of  Eq.  (3)  are  squared  to  give  the  middle  equality  in  Eq.  (7).  This 
equation  is  then  multiplied  by  ^x(/-l,t)^  to  yield  the  lower  equality  in  Eq.  (7).  Each 
equality  in  Eq.  (7)  is  written  in  terms  of  the  parameters  defined  in  Eq.  (5). 
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T-l  T-l  T-l 

'^x{i,t)x{i  -  \,t)  =  ^^a{i,t)x{i-\,t)  +  b{iy^x{i-\,t)^  =r,-  +b{i)zi 

t=0  t=0  t=0 

T-\  T-\  T-\  T-\ 

'^x{i,t)^  =  +2b{i)'^a{i,t)x{i -\,t)  +  b^{i)'^x{i 

t=0  t=0  t=0  t=0  (7) 

=  Si+2biiyi+b^ii)zi 

T-l  T-l  /  N. 

X S (0^/  j  =  +  2/)(z )r,z.  +b'^{i)zy 

t=0  t=0 


Equation  (6)  can  now  be  re-written  as  follows. 


T-l 

^x{i,t)x{i-\,t) 


1 


t=(l 

[tm 

^x(z,0^ 

Z=0  Z=0 


(8) 


To  obtain  the  desired  eorrelation,  the  eorreet  seale  faetor  must  be  ehosen  by 
solving  Eq.  (8)  for  b{i).  To  solve  for  b{i),  both  sides  of  Eq.  (8)  are  squared,  the  resultant 
equation  is  rearranged  into  a  quadratie  polynomial,  and  solved  as  follows. 
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If  the  desired  eorrelation  is  positive  (p(/)>  0),  the  maximum  value  of  h(z)  is  ehosen;  if 
the  desired  eorrelation  is  negative  (/?(/)<  0),  the  minimum  value  of  h(z)  is  ehosen.  In 
the  diseussion  given  below,  the  argument  of  the  radical  is  shown  to  be  non-negative. 

At  this  point,  the  correlation  part  of  the  line  scan  problem  is  solved.  With  h(z) 
calculated  as  given  in  Eq.  (9)  and  following  the  procedure  outlined  above  culminating  in 
Eq.  (3),  simulated  acoustic  noise  signals  can  be  created  with  a  correlation  coefficient 
distribution  which  matches  the  desired  distribution. 

C.  Matching  the  Maximum  Extreme  Value  Distrihution 

A  secondary  point  of  emphasis  was  to  force  the  simulated  noise  to  match  the  MEV 
distribution  for  the  measured  noise.  When  using  a  normal  random  number  generator 
approach  to  create  simulated  acoustic  noise  signals,  the  mean  of  the  MEV  distribution  is 
controlled  by  the  standard  deviation  of  the  random  number  generator.  As  discussed 
below,  the  mean  of  the  MEV  distribution  was  handled  at  the  end  of  the  correlated  noise 
creation  step  by  scaling  each  signal  to  control  the  average  signal  energy.  The  shape  of 
the  MEV  distribution  depends  on  the  correlation  between  points  in  each  A-scan;  in  other 
words,  the  shape  is  controlled  by  the  autocorrelation  function  or  average  power  spectrum 
of  the  noise.  As  described  above,  the  initial  set  of  uncorrelated  noise  signals 
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i  =  \,N)  are  filtered  to  foree  an  average  magnitude  speetrum  whieh  matehes  that 
of  the  measured  noise.  The  summation  proeess  depleted  in  Eq.  (3)  alters  the  average 
speetrum  for  the  output  eorrelated  signals  i  =  \,N).  Unfortunately,  eorreeting  the 

average  speetrum  of  the  eorrelated  signals  baek  to  the  desired  speetrum  would  alters 
signal  eorrelations  and  assoeiated  eorrelation  eoeffieient  distribution;  thus,  motivating  the 
implementation  of  an  iterative  approaeh.  The  goal  is  to  filter  the  noise  going  in  to  the 
eorrelation  proeess  so  that  the  noise  eoming  out  of  the  eorrelation  proeess  will  have  the 
desired  average  power  speetrum  and,  therefore,  the  desired  MEV  distribution  shape. 

Additional  notation  is  needed  to  deseribe  the  iterative  approaeh.  The  output  of  the 

eorrelation  proeess  for  the  iteration  will  be  denoted  i  =  \,N  with 

average  power  speetrum  and  magnitude  speetrum  given  by  |Xy^(/)|^  and  |Xy^(/)|, 
respeetively.  The  filters  will  be  denoted  Fy^(/)  with  the  initial  filter,  Fi(/),  based 
direetly  in  the  magnitude  speetrum  for  the  measured  noise:  F^  (/)  =  |X^  (/ )| .  The  initial 
filter,  Fi(/),  is  the  filter  used  in  ereating  the  uneorrelated  signals  (a  =  a{i,t)  i  =  l,N). 
The  goal  is  then  to  make  the  output  speetrum,  |Xy^(/)| ,  equal  to  the  desired  speetrum, 
\Xin  (/ )| ,  to  within  some  aeeeptable  error.  Defining  to  be  the  sum  of  squared  errors 
between  |X^(/)|  and  |X„(/)|  over  some  frequeney  range  and  to  be  the  aeeeptable 
error  level,  the  iterative  proeess  will  eontinue  until  Ex^^x-  The  filter  is  updated  after 
eaeh  iteration  based  on  the  ratio  of  the  desired  and  output  speetra:  |X^  (/ )|/|X^  (/ )| . 
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For  a  line  scan,  the  iterative  approach  used  to  achieve  the  desired  correlation 
coefficient  distribution  and  frequency  spectrum  can  be  summarized  in  algorithm  form  as 
follows: 


0. 

1. 

2. 

3. 

4. 

5. 

6. 


generate  a  set  of  white  noise  signals  :  i  =  \,N 

apply  the  fdter,  to  the  white  noise  to  yield :  a  =  a{i,t)  i  =  \,N 

calculate  :  xj^{\,t)  =  a{\,t)  xj^{i,t)  =  xj^{i-\,t)  +  b{i)a{i,t)  i  =  2,N 

ff 

calculate  |X^(/)|  and  the  error :  =  ^(]x^(/)|-|X^(/)|f 

fi 


x  =  x^;  stop; 


if  calculate  an  updated  fdter :  k  =  k  +  1 

return  to  step  1 . 


ri(/)  =  Fn(/) 


x„(/) 

Xi-i(/) 


Empirical  evidence  shows  that  the  final  filter  shape  is  strongly  dependent  on  the 
desired  magnitude  spectrum  and  weakly  dependent  on  the  desired  conditional  correlation 
coefficient  distributions.  As  such,  for  a  given  desired  magnitude  spectrum,  we  generally 
run  through  the  iterative  process  one  time  to  establish  the  appropriate  filter  shape.  This 
filter  is  then  used  for  all  additional  runs  to  generate  correlated  noise  signals,  regardless  of 
the  desired  output  correlation  coefficient  distributions. 

Finally,  the  signals  are  scaled  to  achieve  the  desired  mean  MEV,  that  is,  the  mean 
MEV  associated  with  the  measured  acoustic  noise  signals,  denoted  indicated 

earlier,  when  generating  uncorrelated  noise  signals  (a.  =  a{i,t)  i  =  l,N)  using  a  standard 
normal  random  number  generator  with  cr  =  1  will  force  the  mean  MEV  for  the  signals  in 
a  to  unity.  The  mean  MEV  for  the  correlated  noise  signals  ( x  =  x{i,  t)  i  =  \,N)  can  also 
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be  forced  to  1  by  making  the  energy  of  each  correlated  noise  signal  equal  to  the  energy  of 
the  associated  uncorrelated  signal.  The  desired  mean  MEV  can  then  achieved  by  scaling 
by  the  mean  MEV  of  the  measured  noise,  denoted  This  can  be  done  in  the 

computational  loop  or  as  a  post-processing  step  as  follows:  x  =  x{a^laj^)m„^y  where 
cr^  and  cr^  are  the  standard  deviation  of  the  uncorrelated  and  correlated  noise, 

respectively.  Note  that  scaling  the  signals  does  not  change  the  correlation  coefficient 
distribution  since  the  correlation  coefficient  is  scale  independent. 

We  close  the  section  by  showing  that  b{i)  is  purely  real  and  by  considering  some 
special  cases.  To  prove  that  b(i)  will  never  be  imaginary,  the  value  under  the  square  root 

must  be  shown  to  be  non-negative.  There  are  three  values  under  the  square  root  p^(i), 

{\-p^{i)) ,  and  (5,.z.  -r^)  .  Since  p  is  always  between  -1  and  1,  the  first  two  quantities 
are  clearly  greater  than  or  equal  to  zero.  To  show  that  the  third  quantity  is  always 
positive,  consider  the  following  manipulations  of  Eq.  (6?),  noting  that  p^{i)>0  . 

=  ^  ^  n  (10) 

A  few  special  cases  should  be  examined.  The  first  case  is  where  the  desired 
correlation  is  ±  1 .  Erom  Eq.  (9),  it  can  be  seen  that  setting  the  desired  correlation  to  ±  1 
makes  the  denominator  under  the  square  root  equal  to  0,  forcing  b(i)  to  ±oo.  With 
reference  to  Eq.  (3),  this  makes  sense  because  /?(/)  =  ±1  implies  x{i,t)=±x{i-\,t)  which 
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can  only  be  approximated  with  b(i)  =  ±<x)  {a{i,t)j^0).  If  the  desired  correlation  were 
±  1 ,  x{i,  t)  should  simply  be  set  to  ± x(i  -\,t) . 

A  second  special  case  would  arise  if  the  original  eorrelation  (Eq.  (5))  were  equal 
to  ±  1 .  For  this  ease,  rf  =  s^z.  and  the  numerator  under  the  radical  in  Eq.  (9)  goes  to 
zero.  This  eliminates  the  influenee  of  the  desired  correlation  value,  p(i),  and  results  in 
b{i)  =  -rilzi  .  We  note  from  Eq.  (8?)  that  b{i)  =  -ril  Zi  ^  r,-  =-h(/)z,-  will  merely  set  the 
output  correlation  to  0.  However,  when  dealing  with  randomly  generated  noise,  it  is 
exceedingly  unlikely  that  we  will  eome  aeross  two  signals  that  have  a  correlation  of  1 .  If 
this  case  were  to  oceur,  the  problem  could  be  solved  by  generating  a  new 
D,  Two  Dimensional  Generation  of  Spatially  Correlated  Noise 

The  procedure  for  generating  eorrelated  A-seans  which  simulate  an  xy  raster  scan 
follows  directly  from  the  procedure  for  a  line  scan.  Figure  8  defines  some  of  the  notation 
used  in  this  seetion.  For  the  two  dimensional  case,  the  correlated  signals  are  ereated  in  an 
inner  iterative  loop  with  a  second  outer  iterative  loop  used  to  simultaneous  satisfy  the 
eorrelation  distribution  and  frequeney  eontent  requirements.  We  again  begin  by 
assuming  that  a  set  of  data  is  taken  by  measuring  backscattered  signals  at  A  x  M  equally 
spaced  measurement  positions  (see  Fig.  1).  The  correlation  coefficients  between 
neighboring  signals  can  be  used  to  define  the  overall  and  conditional  correlation 
eoefficient  distributions.  These  distributions  ean  be  used  to  generate  N(M-1) 
eorrelation  eoeffieients  between  signals  in  a  row,  denoted  Pf.(i,j),  and  {N -\)M 
eorrelation  eoeffieients  between  signals  in  a  eolumn,  denoted  Pf,{i,j).  These  are  the 
desired  correlation  values  that  will  be  used  in  determining  the  scale  faetors,  similarly 
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denoted  and£>^(/,7).  As  with  the  line  scan  case,  desired  correlation  coefficients 

can  only  be  drawn  from  distributions  conditioned  on  correlation  coefficient  values 
between  previously  generated  signals. 

The  procedure  again  starts  with  the  generation  of  a  set  of  uncorrelated  acoustic 
noise  signals,  now  written  into  a  thee-dimensional  matrix:  2L  =  a{i,j,t)i  =  \,N  j  =  l,M  . 
The  first  row  and  the  first  column  are  treated  as  line  scans.  The  first  output  signal, 
is  set  equal  to  a(l,l,t).  The  procedure  described  above  for  a  line  scan  is  then 
used  to  find  the  scale  factors,  bf.{l,j)  j  =  2,M  and  b^{i,\)  i  =  2,N ,  and  the  associated 
simulated  signals,  j  =  2,M  and  x{i,\,t)  i  =  2,N . 

The  remainder  of  the  correlated  signals  are  created  in  an  iterative  fashion  starting 
with  x(2,2,t).  Each  new  signal  is  simultaneously  forced  toward  the  desired  correlation 
with  two  neighboring  signals.  In  general  notation  (see  Fig.  12),  x{i,j,t)  is  iteratively 
forced  toward  the  desired  correlation  with  x{i-\,j,t)  and  x{i,j -\,t)  by  using  the  sum  of 
appropriately  scaled  versions  of  these  two  signals  plus 

To  generate  the  signal  x{i,j,t),  the  algorithmic  loop  can  be  summarized  as 
follows: 

1 .  Calculate  the  scale  factors  required  to  force  the  desired  correlation  between 
x{i,j,t)  and  x{i-\,j,t),  denoted yC'^(/,7),  and  between  x{i,j,t)  and  x{i,j-\,t), 

denoted  (/,  j) : 


20 


(11) 


2.  Calculate  the  output  signal  using  both  seale  faetors  and  all  three  signals: 


x{i,  j,t)  =  a(i,  j,t)  +  br  {i,  j,  )x(i,  7  - 1,  t)  +  (b  j)x{i  -1,7,0  (12) 


3.  Caleulated  the  aetual  eorrelation  eoeffieient  between  x{i,j,t)  and 
x{i -1,7, t) , denoted  (/, 7) ,  and  between  x{i, 7, t)  and  x(/,  j 
denoted  Pj.  (z,  7) , 


T-l 

^x(z 


X  ^(h  2)  0^  J  X  2)  0^ 


(13) 


T-l 

^x(z  ,7,04h7-l,0 


J^x(i,j,tfl^x(i,j-l,t) 


z=0  V  z=0 
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4.  Calculate  the  eorrelation  coefficient  error  and  compare  with  the  aceeptable  error 
level,  Ep'. 

5pc[Ui)  =  \pc[hi)-pc[Uj\  5pr[ij)=\pr[Ui)-pAUj\  (14) 

If  dp^{i,j)<Sp  and  dpj.{i,j)<Sp  then  stop;  else  go  to  step  5. 

5.  Use  the  current  output  signal  as  a  new  starting  signal:  a{i,  j,  t)  =  x{i,  j,  t).  Return 
to  step  1 . 

The  loop  is  repeated  until  NM  -  {N  +  correlated  acoustic  noise  signals  have 
been  generated.  Note  that  x{i,j,t)  is  correlated  with  eaeh  of  its  four  neighbors  sinee 
x{i,j,t)  is  ealculated  based  on  its  eorrelation  with  x{i-\,j,t)  and  x{i,j-\,t),  and 
x{i  +  \,j,t)  and  x{i,j  +  \,t)  are  ealeulated  based  on  their  eorrelations  with  x{i,j,t)  (see 
Fig.  8). 

Finally,  the  MEV  distribution  shape  and  position  are  addressed.  As  with  the  line 
scan,  an  outer  iterative  loop  is  used  to  force  the  generated  signals  to  have  the  desired 
frequency  content,  and  thus  the  MEV  distribution  shape,  while  maintaining  the  desired 
correlation  coefficient  distribution.  This  iterative  process  follows  direetly  from  the  steps 
outlined  for  a  line  sean  in  the  previous  seetion  with  only  notational  changes  required  to 
aceount  for  three-dimensional  rather  than  two-dimensional  matrices.  The  position  of  the 
MEV  distribution  for  the  simulated  signals  is  again  dealt  with  by  scaling  of  the  simulated 
signals,  following  the  approach  outlined  for  the  line  scan:  x  =  )f^mev  ■ 
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IV.  RESULTS 


A.  Implementation 

The  MEV  distributions  (Fig.  1)  and  conditional  correlation  coefficient 
distributions  associated  with  the  overall  distributions  shown  in  Fig.  3  were  used  as  the 
desired  distributions  to  demonstrate  implementation  of  the  noise  generation  approach.  A 
periodogram  approach  was  used  to  estimate  the  average  magnitude  spectrum  of  the 
measured  noise  for  use  in  creating  acoustic  noise  signals  with  the  desired  MEV 
distribution  shape.  Desired  correlation  coefficients  were  drawn  from  interpolation-based 
conditional  probability  density  functions  which  were  based  on  conditional  histograms 
constructed  from  the  measured  noise.  Spatially  correlated  acoustic  noise  A-scans,  each 
51  points  long,  were  generated  to  simulate  a  50  x  50  raster  scan  at  a  digitization  rate  of 
100  MS/s  (10  ns/point). 

B,  Results 

We  begin  by  showing  an  example  of  the  evolution  of  correlation  coefficients  and 
signals  toward  the  desired  result  for  the  middle  distribution  shown  in  Fig.  3.  The  two 
dashed  lines  in  Fig.  9  show  examples  of  the  how  the  actual  correlation  values, 
and  move  toward  the  desired  values,  Pf.{i,j)  and  Pf,{i,j),  with  each  iteration. 

Figure  10  shows  how  one  uncorrelated  noise  signal,  a{i,j)  (solid  lines,  upper  graphs), 
can  start  at  an  arbitrary  correlation  coefficients  with  its  neighboring  signals  (dashed  lines 
in  the  upper  and  lower  graphs)  and  then  after  iteration  evolve  to  a  correlated  noise  signal, 
x{i,j)  (solid  lines,  lower  graphs),  showing  the  desired  correlations  with  its  neighbors. 

The  method  outlined  for  creating  correlated  noise  was  tested  for  each  of  the 
overall  correlation  coefficient  distributions  shown  in  Fig.  3.  Given  the  stochastic  nature 
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of  this  problem,  each  time  a  set  of  correlated  noise  signals  is  generated,  the  resultant 
histograms  represent  one  realization  of  noise  generation  process.  In  this  section,  we 
show  example  histograms  for  one  realization  of  the  process,  and  we  give  chi-squared 
analysis  results  that  quantify  the  average  agreement  between  the  desired  and  output 
histograms. 

Desired  and  actual  histograms  are  given  in  Fig.  11  for  overall  correlation 
coefficient  distributions  and  in  Fig.  12  for  the  MEV  distributions.  The  histograms  in  each 
figure  show  the  desired  correlation  distributions  (Fig.  3)  and  MEV  distributions  (Fig.  1). 
The  star  data  points  superimposed  on  the  histograms  represent  the  bin  values  for  the 
histograms  determined  from  the  simulated  noise.  For  each  distribution,  signals  were 
generated  to  simulate  a  50  x  50  scan,  yielding  2500  signals,  and  4900  nearest-neighbor 
correlations.  The  success  of  the  approach  is  witnessed  by  the  agreement  between  the 
histograms  for  measured  and  created  noise  signals.  Again  note  that  each  set  of  simulated 
signals  will  give  a  slightly  different  correlation  coefficient  histograms  and  MEV 
distributions. 

The  results  were  also  evaluated  quantitatively  using  a  chi-squared  test.  The  test  is 
based  on  the  agreement  between  histogram  bin  values  associated  with  the  generated  and 
measured  noise,  with  the  measured-noise  based  histograms  representing  the  desired 
values.  Using  the  vector  and  to  hold  the  bin  values  for  the  calculated  and 

measured  noise,  respectively,  the  chi  squared  error  between  the  two  histograms  can  be 
written  as  follows: 


h 


m 


(15) 
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where  k  represents  the  number  of  bins.  This  statistie  was  looked  up  in  a  standard  ehi 
squared  table  with  k-\  degrees  of  freedom.  The  resulting  p-value  is  the  probability  that 
the  same  error  would  result  if  both  histograms  were  indeed  from  the  same  distribution. 
Thus  if  we  piek  a  eutoff  value,  sueh  as  0.05,  any  p  value  above  this  eutoff  is  deemed 
aeeep table. 

This  proeess  was  done  on  10  sets  of  20  x  20  of  simulated  noise.  For  the  MEV 
histogram  shown  in  Fig.  1,  the  average  p-value  was  0.325,  with  a  standard  deviation  of 
0.233.  For  the  overall  eorrelation  histogram,  the  average  p-value  was  0.557,  with  a 
standard  deviation  of  0.340.  For  the  eonditional  eorrelation  histogram,  the  average  p- 
value  was  0.220  with  a  standard  deviation  of  0.404.  These  quantitative  results  indieate 
that  on  the  average  for  MEV,  overall  eorrelation,  and  eonditional  eorrelation  distribution, 
the  distributions  for  the  simulated  noise  are  eonsistent  with  the  desired  distributions  based 
on  the  measured  noise. 

V.  SUMMARY 

An  approaeh  has  been  outlined  for  generating  simulated  aeoustie  noise  with  a 
eorrelation  eoeffieient  distribution  and  MEV  distribution  whieh  matehes  those 
distributions  for  measured  aeoustie  noise.  With  this  approaeh,  a  limited  number  of 
measured  signals  ean  be  used  to  establish  the  eorrelation  eoeffieient  and  MEV 
distributions  whieh  drive  the  eomputer  generation  of  a  large  number  of  simulated 
aeoustie  noise  signals.  Signals  simulated  in  this  manner  are  being  used  in  the 
development  of  a  seale  independent  eorrelation  based  approaeh  to  defeet  deteetion.*^^’^^ 
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1.  List  of  Figures 


Fig,  1.  Histograms  of  maximum  extreme  values  for  1.5  mm  (upper)  and  0.5  mm  (lower) 
measurement  position  spaeing. 

Fig.  2.  Representation  of  measurement  positions  and  assoeiated  data  matrix. 

Fig,  3,  Histograms  of  spatial  eorrelation  eoeffieients  and  assoeiated  fits  for  2  mm 
(upper),  1  mm  (middle),  and  0.5  mm  (lower)  measurement  position  spaeing. 

Fig.  4.  Variations  in  spatial  eorrelation  eoeffieients  along  a  line  sean  for  0.5  mm  (upper) 
and  1  mm  (lower)  measurement  position  spaeing. 

Fig.  5.  Stem  plots  of  overall  and  eonditional  eorrelation  eoeffieient  distributions.  The 
eonditional  distribution  is  for  an  adjaeent  eorrelation  eoeffieient  of  0.58. 

Fig,  6,  Steps  used  in  ereating  spatially  uneorrelated  aeoustie  noise  with  the  desired 
frequeney  eontent:  (a)  white  noise  signal  generated  using  a  standard  normal  random 
number  generator;  (b)  speetrum  of  white  noise  signal;  (e)  filter  determined  based  on  the 
average  magnitude  speetrum  of  the  measured  aeoustie  noise  signals;  (d)  speetrum  after 
filtering;  and  (e)  inverse  Fourier  transform  of  filtered  speetrum. 

Fig,  7.  Representation  of  the  notation  and  steps  involved  in  ereating  a  simulated  line 
sean  of  spatially  eorrelated  signals. 

Fig.  8,  Representation  of  the  notation  and  steps  involved  in  ereating  a  simulated  raster 
sean  of  spatially  eorrelated  signals. 

Fig.  9.  Example  of  the  progression  of  row  (upper)  and  eolumn  (lower)  eorrelation 
eoeffieients  toward  the  desired  values  during  the  iterative  proeess. 

Fig,  10,  Example  of  a  signal,  initially  having  arbitrary  eorrelations  with  adjaeent  signals, 
evolving  during  the  iterative  proeess  to  have  the  desired  eorrelations  with  the  adjaeent 
signals.  The  top  row  shows  the  original  signal  to  be  modified  (solid  lines)  along  the 
adjaeent  signals  (dashed  lines).  After  the  iterative  proeess,  the  final  modified  signal 
(solid  line)  is  shown  with  the  adjaeent  signals  (dashed  lines). 

Fig.  11,  Example  histograms  of  spatial  eorrelation  eoeffieients  for  measured  and 
simulated  noise  signals  at  2  mm  (upper),  1  mm  (middle),  and  0.5  mm  (lower) 
measurement  position  spaeing. 

Fig,  12,  Example  histograms  of  maximum  extreme  values  for  measured  and  simulated 
noise  signals  at  1.5  mm  (upper)  and  0.5  mm  (lower)  measurement  position  spaeing. 
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